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ABSTRACT 

Through extended integrations using the recently-installed deep depletion CCD on the red arm 
of the Keck I Low Resolution Imaging Spectrograph, we present new measurements of the resolved 
spectra of 70 morphologically-selected star-forming galaxies with %ab < 24.1 in the redshift range 
1 < z < 1.7. Using the formalism introduced in Paper I of this series and available HST ACS images, 
we successfully recover rotation curves using the extended emission line distribution of [O II] 3727 A 
to 2.2 times the disk scale radius for a sample of 42 galaxies. Combining these measures with stellar 
masses derived from HST and ground-based near-infrared photometry enables us to construct the 
stellar mass Tully-Fisher relation in the time interval between the well-constructed relation defined at 
z~l in Paper I and the growing body of resolved dynamics probed with integral field unit spectrographs 
at z>2. Remarkably, we find a well-defined Tully-Fisher relation with up to 60% increase in scatter 
and zero-point shift constraint of AM, = 0.02±0.02 dex since z ~ 1.7, compared to the local relation. 
Although our sample is incomplete in terms of either a fixed stellar mass or star formation rate limit, 
we discuss the implications that typical star-forming disk galaxies evolve to arrive on a well-defined 
Tully-Fisher relation within a surprisingly short period of cosmic history. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: kinematics 
and dynamics — galaxies: spiral 



1. INTRODUCTION 

Excellent progress has been made over the past 15 
years in obser vationally charting the history of co smic 
star formation ( |Hopkins fe Beacom|200"6l|Ellis|2008| . At- 
tention has now moved to inter preting the stellar mass 
assembly history of galaxies (e.g., Bundy et al.|2006 ) in a 
framework whereby cold dark matter h alos m erge under 
gravitational instability (White & Rces 1978). Combin- 
ing dynamical and stellar masses tor selected galaxies at 
various look-back times provides insight to the relative 
contributions of dark and baryonic matter through time. 
For the most abundant population of star forming galax- 
ies, the Tully-Fisher relation, which probes the relation- 
ship between luminosity and rotational velocity, is the 
most effective tool in assessing the angular momentum 
dev elopment of disks as they assemble in dark matter ha- 
los dFall fc Efstathiou||l980l iBlumenthal et al.Hi984| [Mo 
et al.| | 1998[ |lNavarro L Steinmetz||2000| ). Of particular 
value is the stellar mass Tully-Fisher relation (M„-TFR), 
which is tighter than a ny scaling relation co nfined to a 
given luminosity band ( |Bell fe de Jong||2001| making it 
highly sensitive to self-sim ilar dis k growth through time. 

In Paper I of this series ( |Miller et alJ|2rJi~l~l ) we demon- 
strated, using extended integrating times with DEIMOS 
on the Keck 2 telescope together with an optimal tech- 
nique for extracting kinematic rotation curves to a loca- 
tion corresponding to 2.2 disk scale radii (2.2r s or ^.2), 
a stellar mass TFR at z ~ 1 that has virtually the same 



intrinsic scatter as that observed locally. This indicates 
that most disk galaxies at LMC-mass or higher are dy- 
namically mature by z ~ 1 , in contrast to the conclusions 
of earlier studies based on sho rter exposure times or other 



velocity extraction methods dConselice et al. 2005 Flo 



res et al.|2006| Kassin et al. |2007 ). Moreover, the zero 
point shift in the relationship from z ~ 1 to z ~ 0.2 is 
marginal (AAf* ~ 0.04 ± 0.07 dex), which we found was 
due to an evolving mix of the baryonic components that 
likely dominate within r 2 .2- While the exact fraction of 
dynamically mature disk galaxies at a given redshift is 
debated, the results of Paper I for disks at z ~ 1 con- 
trasts markedly with the dispersion-dominated dynamics 
of massive, s tar-forming galaxies observed 2-3 Gyr ear- 
lier at z > 2 |G cnzcl et al. 2006 ; Forster Schreiber et al 



2009| |Jon es et al. 2010). In a ddition to a large increase 
Cresci et al.|(| 2009| ) show the z ~ 2.2 M„-TFR 
0.41 ± 0.11 dex from that found 



m scatter, 
has an offset of AiW« 
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in Paper I. The implication is that star-forming galax- 
ies very rapidly establish their present dynamical state 
during the intervening period. Thus, resolved dynam- 
ical data in the redshift interval 1 < z < 2 is needed 
to make progress, and large samples will be valuable to 
understand the demographic trends. 

There are several observational challenges in pursuing 
the M^-TFR beyond z ~ 1. Multi-object capabilities 
are essential in creating a large sample to address ques- 
tions such as those posed above, but these capabilities 
are only currently available at optical wavelengths. The 
emission line of choice, [O II] 3727 A moves to a wave- 
length region of low detector sensitivity, and the physi- 
cal scale being probed is comparable to just a few times 
the average seeing (a diameter of 5r s at z ~ 1.6 corre- 
sponds to ~2.5 arcsec). This second paper in our series 
was motivated by the installation of a deep depletion, 
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red-sensitive CCD on the red arm of the Low Resolution 
Imaging Spectrograph (LRIS, Oke et al. 1995[ ) which has 
significantly improved the efficiency in the wavelength in- 
terval 8,000 A to 10,000 A corresponding to [O II] in the 
redshift interval 1.0 < z < 1.7 (Fig. [T]). This allows us to 
extend our earlier work, to try to understand the tran- 
sition from the dispersion-dominated systems observed 
above z~2, and to evaluate the nature of the M*-TFR 
2 Gyr earlier than that probed in Paper I. We present 
the results based on a sample of 70 morphologically- 
selected star-forming galaxies for which we have obtained 
the necessary extended integration times, following the 
techniques developed in Paper I. 

The plan of the paper is as follows: in §2, we describe 
the sample selection criteria, the LRIS data and the HST 
Advanced Camera for Surveys (ACS) resolved photome- 
try and stellar mass estimates; §3 explains our technique 
in determining the fiducial rotation velocities, as well as 
various tests of the chosen method; in §4, we present our 
results of the TF relation from 1.0 < z < 1.7, and in §5 
we discuss our results in light of current disk assembly 
theory and pioneering observations to even higher red- 
shift. Throughout the paper we adopt a J7a = 0.7, £l m 
— 0.3, Ho = 70 km sec -1 Mpc -1 cosmology. All magni- 
tudes refer to the AB system. 

2. DATA 

We selected our sample of spectroscopic targets in a 
similar manner to that in Paper I so as to span a broad 
mass and luminosity range. A more detailed look at sam- 
ple representativeness can be found in the Appendix, 
but the salient points are summarized in the following 
section. The necessity of HST ACS imaging for accu- 
rate disk size measurements led us to focus on five fields 
with multi-wavelength coverage for selecting targets: (1) 
the Great Observatories Origins Deep Survey (GOODS) 
North and (2) GOODS South fields ( [Giavalisco eTaf] 
20041, (3) the Small Selected Area 22 (SSA22) field SHM 
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et^||1991||Chapman et al.||2004| [Abraham et al^2UU7|), 
(4) the Extended Groth Strip (EGS) field (jDavis et al. 



Strip (EGS) 



IP 



2005)) , and (5) the Cosmic Evolution Survey (COSMOS) 

A summary of masks 



field ( Scoville et al. 



gets, spectroscopic o 



2007). 



tar- 

jservations and available space and 
ground-based imaging in each field can be found in Table 

ED 

In selecting our targets we applied an apparent mag- 
nitude limit of %ab < 25 and selected sources in the red- 
shift range 1.0 < z < 1.7 (using photometric redshifts 
when a spectroscopic redshift was unavailable). This led 
to a sample of 50 targets per LRIS pointing. Poten- 
tial targets were then visually inspected by at least two 
of us (RSE, KB, SHM) on the ACS images and poten- 
tial targets with compact morphology, without diffuse 
extension, were given low priority but not excluded^] 
From a prioritized list of typically 25 suitable targets per 
LRIS pointing, on average ~10 of these were included per 
mask due to the requirements of both suitably-oriented 
slits and simultaneous accommodation of unrelated tar- 

5 In considering the effect of excluding unresolved or compact 
sources, it is important to recognize that the dynamical measure- 
ments for spatially unresolved galaxies will likely add to the ob- 
served scatter through their measurement uncertainties rather than 
because those sources have an intrinsically larger scatter around the 
relationship. 
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Fig. 1. — The redshift distribution of 70 galaxies in our new 
sample compared to that in Paper I. Line emission was detected in 
most of the z > 1 galaxies targeted as shown in solid black. The 
dotted orange line shows the distribution for those galaxies with 
no detected emission. 

gets for other science campaigns. In designing multi-slit 
masks for LRIS, we aligned 1 " slitlets to the major axes 
of ou r targets as fit by SEXTRACTOR (Bertin & Arnouts 
19961. The position angle (PA) of the mask was then 
selected to minimize slitlet tilts while maximizing target 
count, and we always ensure that target slitlets are within 
±30° of the mask PA. We used the 600 1 mm -1 grating 
blazed at 10000 A, providing a instrumental velocity res- 
olution of 58 km s _1 at 9000 A thereby exploiting the 
sensitive region of the newly-installed CCD. The spectro- 
scopic LRIS data were reduced b y one of us (AN) using 
techniques described in det ail by Newman et al. (2010) 
with the code developed by Kelson (2003]) , in which sky 
modeling and subtraction is carried out with sub-pixel 
sampling of the background. The seeing, which varied 
from 0'.'78 to l'/05, was measured by taking the weighted 
average of the FWHM of the alignment stars on each 
mask. 



2.1. Photometric Components 

With the existing HST ACS data (see Table[l]for avail- 
ability per field), we conduct our measurements of disk 
radii and inclination of each galaxy using the reddest 
available filter (F814W or F850LP). We ensure that light 
from the bulge component does not confuse the fitted 
disk component by conducting a bulge-disk decomposi- 
tion, where appropriate, similar to our procedure in Pa- 
per I. 

The HST images were anal yzed for pho tometric com- 
ponents using the GALFIT3 (Peng 2010) least squares 
elliptical-fitting method described in detail in Paper I to 
which the reader is referred. Briefly, we fit an exponen- 
tial disk component plus a de Vaucouleurs' bulge profile 
to every galaxy. Those galaxies which yielded unphysical 
solutions were re-fit with a single Sersic profile compo- 
nent, where the Sersic index (n) was allowed to vary. 
Such cases generally represent disk galaxies which are 
bulgeless and/or more clumpy and irregular than regu- 
lar well-formed spirals. For the present sample, ~63% of 
our galaxies were fit using a 1-component n-varying Ser- 
sic profile fit, and ~37% were adequately fit with a two- 
component bulge and disk solution. This mix provides 
a valuable indication of the morphological distribution 
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TABLE 1 

Summary of LRIS observations and data 



N obj b 



DATES OF OBS 



HST Filters 



egs 2 11(3,3) 30.0 ks 26-28 ju„ 2009 1'.'05 

12(2,2) 14.4 ks " tf.'93 

SSA22 2 13(3,3) 22.5 ks " CK'78 

12(2,1) 12.0 ks " / /89 

GOODS N 1 8(1,3) 34.8 ks 5-6 Apr 2010 0'.'86 

COSMOS 1 9(2,1) 28.8 ks 7 & 10 Jam 2011 CK'82 

goods s 1 5(1,1) 25.2 ks 28 Fob 2011 tf'98 



F606W, F814W 
F814W 



Palomar (JK S ) 
UH 2.2m (JHK a ) 



F435W, F606W moircs (K s ) 

F775W, F850LP on Subaru 
F814W 

F435W, F606W 



WIRCam (K s ) 

on CFHT 
ISAAC (K s ) 



F775W, F850LP on bso vlt 



Bundy et aT| | |2006| 



Capak et al| ( |2004} 



Bundy et aT] | |2009| 



McCracken et aI.||20To) 



Retzlaff et aL| ( |2010} 



a number of masks per field 

b number per mask and those of which are (compact, passive), see [3.1 
c total integrated exposure time per mask 
d dates observations were taken at Keck I 
filters of HST ACS imaging available for size measures 

f Ground-based infra-red photometry available for SED fitting, in addition to further optical filters (not shown) 
g Reference to photometry matching catalog 




1 1.0 1.2 1.4 1.6 2 4 6 8 
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z 




02468 10 12 14 2 4 6 8 

r s (kpc) Emission Extent (r/r s ) 

Fig. 2. — Distributions of redshift, apparent i magnitude, stellar mass, exponential disk scale radius (r B ), and the extent of [O II] line 
emission for our sample (measured in terms of the disk scale radius r s ). The sample is partitioned into the three spectral categories 
discussed in the text: extended line emission (black line), compact line emission (green dashed line), which does not extend past the seeing 
or 1 scale radius, and passive sources (orange dotted line). The fiducial radius of 2.2 times the disk scale radius r s is plotted with the 
vertical black dashed line. See [j3]for further details. 
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of our sample, suggesting less than half are well-formed 
disks. For comparison, this mix was 60:40 in Paper I. 

We ran GALFIT using HST data of the reddest avail- 
able filters (noted in Table [l]). Where imaging from more 
than one filter is available, we run GALFIT on each and see 
<10% difference between bands. In our comparisons be- 
tween the F435W filter in Paper I and the redder filters, 
we see no systematic offset, but typically a <5% increase 
in the scatter of the radius fit. We assume an additional 
uncertainty of this nature exists in the measurements of 
this dataset since we are sampling rest-frame blue-UV 
light in the majority of the sample. To achieve conver- 
gence on the GALFIT parameters and to assess their er- 
rors, we run a similar Monte Carlo analysis to that in 
Paper I. We found the parameter output distributions 
were much narrower than the input distributions, sug- 
gesting convergence. Final parameter uncertainties from 
the Monte Carlo distributions are better than 10% on 
average, and the additional uncertainties arc added in 
quadrature to the observational error and formal fitting 
errors. 



2.2. Stell 



ar masses 



Our method of estimating stellar mass in our galax- 
ies fol lows the work initially presented in Bundy ct al. 
( 20051 ), followed by the analysis presented in |Bundy et al. 
(2009). Further details can be found in those papers. 

Briefly, stellar masses are derived using a matched cat- 
alog of multi-band ACS, available optical and ground- 
based infra-red photometry (see Table [I] for ACS and 
infra-red availability per field). A Bayesian code fits 
the spectral energy distribution (SED) derived from 
2" ACS, and ground based optical and infra-red photom- 
etry, adopting the best spect roscopic redshift and this 
SED is compared to a grid of Bruzual & Chariot ( 2003 ) 
models that span a range of nietallicities, star formation 
histories, ages and dust content. The stellar mass is cal- 
culat ed from th e deriv ed -fT-band mass/light ratio assum- 



ing a Chabrier (2003J initial mass function (IMF). The 



probability for each fit is marginalized over the grid of 
models giving a stellar mass posterior distribution func- 
tion, the median of which is the catalogued value. At 
the magnitudes probed in this survey, the uncertainties 
inferred from the median 68% equivalent width of these 
posterior functions is 0.174 dex. 

In order to determine the systematic uncertainties aris- 
ing from combining different sets of photometry into a 
single sample, we compare mass estimates for our sam- 
ple in SSA22 using different combinations of the avail- 
able photometry, e.g. the filters available in EGS. The 
resulting dispersion in the distribution of stellar mass es- 
timates ranges per galaxy from 0.005 to 0.1 in dex, with 
a median standard deviation of 0.052, which we add in 
quadrature to each galaxy's SED median posterior width 
as an additional systematic uncertainty. The combined 
stellar mass uncertainty is better than 0.2 dex in over 
80% of our samplej^] As in Paper I we use the stellar 
mass enclosed within 2.2r s , as scaled from the reddest 
available HST broad-band flux, in order to compare like 

6 This does not include systematic uncertainties regarding our 
incomplete knowledge of stellar populations in the various types 
of galaxies present in this sample at these redshifts — uncertainty 
which affects all present studies of this nature. 



velocity with like mass. The average reduction in stellar 
mass at r2.2 is -0.187 dex compared to the total stellar 
mass estimates to the Kron radius (Kron 1980). For the 
present sample, we successfully probe down to a mini- 
mum total stellar mass of log M„/M Q of 8.73 dex, with 
a mean of 10.04 dex. 

3. DYNAMICAL ANALYSIS 

Modeling rotation curves beyond z ~ 1 is a chal- 
lenging endeavor in many respects. In addition to the 
shift of the diagnostic [O II] 3727 A emission line into 
the sky-dominated red spectral region, the more dis- 
tant galaxies are smaller in angular scale making re- 
solved measurements challenging from the ground with- 
out the use of adaptive optics (AO). Prior to the avail- 
ability of near infra-red (NIR) multi-IFU instruments 
(e.g., KMOS on the VLT), multi-slit spectrographs re- 
tain several advantages for delivering kinematics for large 
samples of galaxies over a wide redshift range. Only 
a fraction of the presently-available IFU datasets have 
been observed with AO, and the typical wider wave- 
length ranges and multi-plexing capabilities of spectro- 
graphs such as DEIMOS, LRIS, and the soon-to-be- 
commissioned MOSFIRE on Keck enable extended ex- 
posures on survey sources without a penalty in spectral 
or spatial resolution. Improved multiplexing of IFUs will 
mitigate this current deficiency as 2-D spectroscopy offers 
clear advantages in optimal fitting of complex velocity 
fields (|Puech et al.|[200 8j |Jones et "afl2010l |Genzel et al. 
2011[[Forster Schreiber et al.||2011[ |Contim et al.||201ip . 
For now, multi-slit spectrographs retain a key advantage 
in survey efficiency. Here we exploit the improved red- 
sensitivity of a newly-installed LRIS deep depleti on CCD 
from Lawre nce Berkeley National Laboratories (Rockosi 
et al. 2010 ) to secure resolved dynamics for 70 galaxies 
in the redshift interval 1 < z < 1.7 with a lower average 
stellar mass range than the many extant IFU samples. 

We discuss our motivation in Paper I to determine the 
rotational velocity at a fixed fiducial radius, and why 
we choose this to be T2.2, or 2.2 times the scale radius 
(2.2 x r s ) of the galaxy, as measured from the broad 
band HST optical filter. A maximum measured velocity, 
V max , is highly dependent on the inconsistent distribu- 
tion of emitting gas in a given galaxy for studies based 
on nebular emission lines such as ours. So it is impor- 
tant instead to base a measurement for the TF relation 
at a consistent location in the disk galaxy such as V2.2, 
the velocity found at T2.2 (further justified in Paper I, 
and references therein). In the following we will briefly 
review our methodology for extracting rotation curves 
and determining the velocity, V2.2 at this fiducial radius, 
as well as testing the uncertainties given the increased 
distance of our sample c.f. the Paper I sample. 

3.1. Rotation Curve Modeling 

As in Paper I we adopt the empirically-motivated arc- 
tan function as the starting point of our modeling pro- 
cedure, viz: 

V = V + -V a arctan( ) , (1) 
7r r t 

where Vb is the central or systematic velocity, ro is the 
dynamic centre, V a is the asymptotic velocity, and 
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Fig. 3. — Examples of our kinematic analysis: (from left to right) the HST image of the galaxy in ACS F850LP or F814W filter; the 
2-D [Oil] da ta in the reduced LRIS spectra; the modeled emission trace using an arctan function blurred by both the measured seeing and 
dispersion (j]3,l[l; the residual obtained by comparing the latter two signals. 
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is the turnover radius, which is the transitional point 
betwe en the rising and flattening p art of the rotation 
curve ( |Courteau||1997[ |Willick|[T999 ) . 

Paramount to modeling emission lines for rotation 
curves is that the emission being modeled is adequately 
resolved. There are two situations in which this may not 
occur: (1) the angular extent of the galaxy is insufficient 
compared to the seeing and spatial sampling of the in- 
strument, and (2) where the galaxy may be sufficiently 
extended in broad band imaging but the spectroscopic 
emission remains unresolved. Our morphological selec- 
tion has effectively de-selected the first category but we 
must now deal with the second. 

Ideally, the line emission must extend beyond the see- 
ing disk. To robustly model the dynamics with the arctan 
function we have generally found that this corresponds 
to a simple criterion that the line emission must extend 
to at least one disk scale radius (r s ), as measured by 
the broad-band photometry. By testing our modeling 
code with simulated data we found in Paper I that too 
many possible models are consistent with the data if this 
simple criterion regarding the emission line extent is not 
met. By progressively removing flux from the outer disk 
emission in the current LRIS data, we have confirmed 
this result from Paper I. We thus applied this criterion 
to define that subset of our sample for which emission 
line data can be appropriately modeled. Studies which 
do attempt to dynamically model unresolved kinemat- 
ics will very likely introduce scatter into their scaling 
relations. Unresolved rotational velocity gradients, in 
both IFU and multi-slit datasets, will propagate through 
the instrument as increased dispersion, degenerate with 
the intrinsic velocity dispersion in the optical and NIR 
-emitting gas of the galaxy. Especially if the unresolved 
velocity gradient of the inner solid body rotation is the 
only part of the galaxy with sufficient signal-to-noise, fit- 
ting a model such as the arctan function in Eqjl] will re- 
sult in poorly constrained extrapolated velocities toward 
the outer disk, e.g. any velocity wishing to sample the 
flat part of the rotation curye for a Tully-Fisher study 
dCourteau) | l997l |Miller et^[2^lT| ). 

(Jt the 71) spectra discussed earlier, 60% (42) show 
emission which is sufficiently extended according to the 
above criterion from which a rotation curve can be mod- 
eled. Of the remainder, 20% (14) have detectable line 
emission but it is too compact, i.e. does not extend above 
r s , and 20% (14) appear to be passive without emission. 
This 60/20/20 composition of extended, compact, and 
passive spectral emission respectively, is similar to the 
equivalent mix (55/25/20) seen in Paper I, and as before, 
these three sub-samples are not statistically different in 
redshift (Fig. [2|. For the spectrally compact sample, 
it appears that 6 out of the 14 galaxies are simply too 
small to resolve, and the other 8 are well extended in 
their broad-band optical disk, but are only undergoing 
significant star formation in their central regions, within 
1 disk scale radius. These spectrally compact objects as 
well as the spectrally passive objects have larger stellar 
masses by 0.2 dex than disks with extended emission. 
Spectrally passive disks have only slightly redder I — K 
colors (by > la) than the distribution width of extended 
emission galaxies, but show no significant difference in 
size compared to those with extended emission. 

As in Paper I (see for full detailed modeling procedure) , 
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Fig. 4. — Three galaxies at z ~ 1 observed with DEIMOS, com- 
pared to a simulated observation at the different spatial and spec- 
tral resolution of the LRIS instrument directly below. From left 
to right, the columns show (i) the spectral data, (ii) the best-fit 
model, and (iii) the residual of the data from the model. 



we blur the arctan function in both the spatial and spec- 
tral dimensions to account for seeing and instrumental 
dispersion, respectively; these are highly degenerate in 
many cases. In order to attempt to break this degener- 
acy we measure the position-dependent dispersion from 
the LRIS spectra in each galaxy, and separately measure 
the seeing from the alignment stars in each slit mask. In 
the model, we implement the dispersion measured from 
the spectrum, reduced by a factor fit to remove addi- 
tional dispersion of the line from blurring by seeing. We 
then blur the model by the seeing, and re-normalize the 
flux across each spatial bin (see Fig (3l for examples). The 
factor to account for the effect of the seeing in the dis- 
persion is fit by the Levenberg-Marquardt least-squares 
solver, along with the parameters of the underlying arc- 
tan function shown in Eq. [I] We then correct the ex- 
tracted values of V2.2 at r 2 .2 for the effects of disk in- 
clination with the minor-to-major axis ratio b/a derived 
from our GALFIT results. 

For 35 galaxies (half our original sample) our spectra 
are sufficiently deep that we can robustly follow the line 
emission to r 2 .2 (after accounting for seeing). For the 
remaining 7 galaxies with extended line emission, we are 
forced to infer V2.2 by extrapolating beyond the point of 
maximum emission extent in the best fit rotation curve. 
The error budget for this modeling method is discussed 
in Paper I, and we discuss further uncertainties regarding 
this specific application of our spectral models below. 
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1.6 1.8 2.0 2.2 2.4 2.6 
DEIMOS log[V 22 /km s" 1 ] 

Fig. 5. — Testing the method: For a subset of galaxies from Paper 
I in the redshift interval 0.85 < z < 1.3 we re-sample their original 
DEIMOS spectra to the spatial and instrument resolution of LRIS 
and compare the derived V2.2 with the original measures (see text 
for details). 

3.2. Testing the Model 

At a redshift z ~ 1.5, a typical scale radius of 4 kpc is 
only ^0.5 arcsec and so clearly seeing and instrumental 
dispersion blur much of the detail in our LRIS spectra. 
As it is not always clear visually that there is an intrinsic 
rotation curve with a characteristic turnover point, we 
need to demonstrate the reliability and precision of our 
modeling technique. 

To facilitate this, we selected galaxies in the redshift 
range 0.85 < z < 1.3 studied with DEIMOS in Paper 
I and re-sampled these data to the spatial and spectral 
resolution of the LRIS data. With the cosmology as- 
sumed, the average angular scale of the LRIS galaxies 
(at (z)=1.31) is 8% smaller than that for our chosen 
DEIMOS sample (at (z)=0.98). To account for differ- 
ences between Paper I and the present survey, we there- 
fore resample for the different spatial pixel scales (0'.'1185 
with DEIMOS to fX'27 with LRIS) and match the LRIS 
spectral dispersions (58 km s _1 at 9000 A with 27.9 km 



lution, the scatter in the resulting estimated V2.2 was as 
large as 1 dex in log[V/km s _1 ] compared to that of the 
arctan-based fits at DEIMOS resolution. We also tested 
whether we could recover V2.2 by measuring dV/dr at in- 
tervals of 0.1 scale radii along the disks, but the scatter 
in the estimated V2.2 was on average 0.5 dex, depending 
on the radial extent of the emission. 

4. RESULTS 

The primary result of this paper, the stellar mass 
Tully-Fisher relation (M*-TFR) for the redshift inter- 
val 1 < z < 1.7, is shown in Figj6j As in Paper I, we 
plot the enclosed stellar mass within r 2 .2 against the de- 
projected rotational velocity V^- To interpret the rela- 
tionship with respect to our earlier work, we fit a linear 
regression using a least-squares method incorporating the 
intrinsic scatter (o~i n t) added in quadrature to that of the 
velocity and accounting for errors in both the ordinate 
and abscissa. As in Paper I, we focus our analysis on 
the inverse fit, assuming velocity as a function of mass. 
However, following traditional convention, the relation is 
displayed with V2.2 on the x-axis and M„ on the y-axis. 

Our linear fit can be written as: 



log 



Mr 



b log 



V 2 . 



km 



log 



M 

M Q 



(2) 



where M = 10 lu Af Q . With a slope b =3.869 fixed from 
that derived for the total sample in Paper I , we find 
that a=1.692±0.060, cr int =0.080 dex in log[V/km a' 1 ], 
and the total rms is 0.117 dex, which is reasonable con- 
sidering the median error is 0.102 dex. In stellar mass, 
this corresponds to <7j„ t =0.310 dex in log[M*/M ] and 
an rms of 0.452 dex (where the median error is 0.186 
dex). 

While we recover a zero-point for the fixed-slope M„- 
TFR similar to that found in Paper I for 0.2 < z < 1.3, 
the scatter in the relation has apparently increased as 
much as 60% as compared to the local relation. Part 
of this increased scatter could be attributed to the more 
challenging observations undertaken with LRIS. We note, 
for example, that the additional scatter introduced by 
using LRIS c.f. DEIMOS as described earlier (~0.098 
dex in log[V/km s^ 1 ]), implies our intrinsic scatter could 
be as low as 0.063 dex. Given the uncertainty in the 



s for each pixel). We then analyzed the resampled scatter is about 0.01 dex, this value is not signi ficant! 



DEIMOS data using the same rotation curve fitting pro- 
cedure as we use for the LRIS analysis, and the results 
can be seen in Figs. [4] & [5j Encouragingly, when we 
subtract the recoverea V2.2 of the resampled DEIMOS 
data from that of the original data, we find a weighted 
mean of less than 0.002±0.009 dex, and the scatter in 
the relation between the two is 0.098 dex, which is simi- 
lar to the rms scatter found in the Tully-Fisher relation 
established locally and in Paper I. Re-sampling objects 
at z ~ 1.0 to the angular diameter distance of that found 
at z <~ 1.3 did not make significant changes in the recov- 
ered V2.2J given this corresponds to only a 8% reduction 
in the angular scale. 

We also tested the data with models simpler than the 
arctan-based model. This addresses a possible concern 
that we might be over-fitting the data given the apparent 
lack of resolved detail at this somewhat higher redshift. 
By assuming a linear fit to the emission at LRIS reso- 



different th an that found in Paper I (0.058 dex, |Miller| 



et al. 2011). So while it is possible that the relation is 
less well-defined by a broadening of up to 60% compared 
toz~ 0, where average scatter is found to b e ^0.05 dex 



(i.e., Pizagno et al.||2005 : |Reyes et al.||2011 |, we cannot 
rule out the possibility that the tightness 01 the relation 
in this study is the same as that found locally. 

We also show in Fig. [6] the stellar masses for two im- 
portant subsets of the originally-targeted sample whose 
resolved dynamics cannot be determined. The first cate- 
gory are the 14 spectrally passive objects which are clearly 
some of the most massive galaxies in the sample. Con- 
ceivably these passive galaxies are evolving off the TF 
relation and transitioning onto the red sequence. On the 
other hand, the stellar masses of the spectrally compact 
objects do not significantly differ from that for main sam- 
ple displaying suitably extended emission. Most of these 
galaxies appear to have extended disks in the HST ACS 
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Fig. 6.— The stellar mass Tully-Fisher relation (M*-TFR) de- 
rived for 42 suitable galaxies in the redshift interval Kz<1.7. The 
least squares inverse fit relation is shown as a solid black line and 
the rms scatter is displayed as a pair of dotted lines. The best-fit 
M*-TFR found from z ~ 0.2 to z ~ 1 in Paper I is shown as the 
dashed grey line (which is almost coincident with the results of this 
study). The rms scatter around the Paper I fit is shown the pair 
of grey dotted lines. Although the zero point has not changed in 
the intervening 2 Gyr, the apparent scatter has increased. We also 
show with arrows the stellar masses of two categories of galaxies 
for which we cannot extract measurements of V2.2 (see §3.1 for 
details). 

images, although star formation appears to have largely 
shut off in their outer disks as evidenced by their spec- 
trally compact nature. Conceivably, these galaxies could 
comprise systems with actively forming bulges, up to a 
few times 10 10 M© in stellar mass. 

Finally, to place our Af*-TFR in the context of ear- 
lier work, Fig. [7] summarizes recent results for the M* 
Tully-Fisher relation over the full redshift range 0<z<3. 
Notwithstanding possible selection biases discussed in §5 
it is illustrative to consider this figure. Evolution in the 
M*-TFR zero-poi nt can be normalize d with respect to 
the local study of Reyes et al. ( 2011| , calibrated from 
the parent sample ot which came from SDSS and who 
also adopt the same V2.2 parameterization used in this 
series. To ens ure an adequate com parison, we re-fit our 
data using the Reyes et al. ( |2011 ) slope, (which we note 
is not significantly differe nt from the slope we find). Be- 
cause [Reyes~et7j^ (2011 ) use the total mass rather than 
that enclosed at r 2 .2, we move their zero-point by -0.187 
dex in stellar mass to account for the difference in the 
enclosed mass as compared to the total stellar mass. We 
make a further co rrection to accoun t for the different 
IMF assumed, since|Reyes et al. (20111) adopted a Kroupa 



IMF, which both sets of authors agree ac counts for 
a 0.05 dex decrease compared to our assumed |Chabrier] 
( 2003|). In total, this is a -0.137 dex calibration between 
tne Reyes et al. (2011) study and our work, resulting in 
only a marginal change in the M„-TFR zero-point since 
z ~ 1.7. Formally, the linear regression shown in Fig. [7] 



yields a stellar mass zero-point change of only 0.02±0.02 
dex over the last 9.8 billion years. 
We also include results at z ~ 2.2 and z ~ 3 from 



Cresci et al. 



(2011) respec- 



(2009) and Gnerucci et al. 
tively, which suggest a surprisingly rapid and significant 
evolution b eyond z ~ 1.7. How ever, the average stellar 
mass in the Cresci et al. (20091 sample is 6.3xl0 10 Mq, 



and range from 0.62-31.6 x 10 Mq, whereas our sample 
average is 1.3xl0 10 M , and ranges from 0.05-17.2 xlO 10 
M Q (both using a Chabrier IMF). 

As a theoretical comparison we present preliminary 
M*-TFR evo lution results computed using Galacticus 
vO.9.1 r631p] ( penson|20 12 ) , which includes various mod- 
ules for star formation and AGN feedback as well as 
state-of-the-art theoretical predictions for baryon-halo 
evolution. Shown in Fig. [7] are Galacticus models for 
both disk-dominated galaxies (bulge-to-total mass frac- 
tion < 0.3) and all galaxies. In each case, the mean zero- 
points correspond to galaxies with stellar masses between 
0.5 and 2xl0 10 Mq. Interestingly, the results from mod- 
els of all galaxies shows a similar, non-evolving trend in 
the M*-TFR zero-point, whereas the the disk-only sub- 
sample of the models depart in evolution to the slow side 
of the observed relation as redshift increases, but never 
more than 0.1 dex — even to z^3. In this figure we 
also in clude the semi-analytical model of |Dutton et al.| 
(20111, which is consistent w ith a sample s ubset of the 
results originally presented in Kassin et al. ( |2007 N=73 
from 544 galaxies). Their projected evolution departs 
from the results at z ~ 1 in this study and Paper I. 

5. DISCUSSION 

The major goal of this paper was to extend our ear- 
lier work in Paper I to higher redshifts to determine the 
rate at which disk galaxies arrive on the Tully-Fisher re- 
lation. Figure [7] presents a puzzling discontinuity in be- 
havior over a relatively narrow period of cosmic history. 
A key uncertainty in the interpretation of this figure, 
however, is that not all of the sample satisfies our cri- 
terion of providing a reliable estimate of V2.2- Although 
we have maintained continuity in our selection techniques 
from Paper I to that in this paper in the redshift inter- 
val z~ 1.0- 1.7, similar criteria of morphological appear- 
ance or extended emission lines are difficult to usefully 
apply to the star-forming galaxies being studied beyond 
z ~ 2. As redshift increases in Fig. [7J some subset of 
a mass-selected or a star formation rate (SFR)-limitcd 
population is likely not included in the TF analysis. 

Such a selection bias may be more or less important de- 
pending on the physical model we wish to test. If disks 
grew from chaotic beginni ngs, where clumps collide at 



the centers of dark halos dvan den Bergh et al.| 1996 



Abraham et al. 1997 Fllmegreen &: Klmegreen "[2005 ) 



then the TF relation will continue to evolve as the chaotic 



behavior and mergers (Hammer et al. 2005) transitions 
into more regular accretion onto well-tormed disks. In 
this scenario, we could interpret the lack of observed evo- 
lution as the result of a bias towards "mature" systems 
at all redshifts. However, since we expect high SFRs in 
this violent phase, these less mature systems are more 

7 computed with an input parameter file which is available from 
http://sites.google.com/site/galacticusmodel/downloads/paramctcr- 
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The left pane l shows the progress in constraining the evolution o f the M«-TFR zero-p oint JConselice et al.| 
2008 jPutton et al.|20li] l before the present series of papers, with exception of Vergani ct al. (2012), which : 
servic e while this paper was under review at ApJ. We note that all of the IFU studies IjFuech et al.|2008| 



120051 |Puech et al -l 

i appeared on the arXiv pre-print 
JVergani et al.|2012||Cresci et al.| 

|2009[) lie suggestively offset from the local relation at the same interval in dex, with exception ol tne mgnest redsnut study done at present 
QUnerucci et ah||2011} . On the right we show the im provement of our r esults, showing both the redshift dependent M*-TFR scatter and 



at z~0 (orange circle 



We compare the res ults from this study 
We also compare 
plot, please see ijiTand 



zero-point, assuming a fixed slope as determined by |Reyes et al.| l |2011| 

(red circle) and Paper I results ( green circles) to tho se a t higher redshift ( jCresd et al.)|2009| [Gnerucci et al ||2011| l 
semi-analytic models (SAM) from |Dutton et ah] (2011) and Benson] (GALACTICUS, |2U12| For a lull discussion ot this pic 

likely to be included in our analysis, unless dust obscu- 
ration or possible nuclear concentration of star formation 
makes extended emission lines difficult to detect (i.e., our 
spectrally compact and passive targets). 

Alternatively if disks form predominantly via smooth 
accretion over time, they could easily evolve along self- 
similar dynamical relations, making selection effects 
much less important in understanding the lack of evo- 
lution. How then could we interpret the more drastic 
offsets impli ed by Fig. FT]? Even wit h smooth, cold ac- 
cretion (e.g., Dekel fc BiFnboim|2006[ ), if galaxies formed 
stars more vigorous ly and erratically from z~2-4 than at 
present tim es (e.g 



tight lock-step. Thus the emergence of a recognizable 
M*-TFR at z~1.5 is already a result of significance - 
even if it only applies to some subset of the star-forming 

population. 

At z^3, Gnerucci et al. (2011) argue that the nascent 
TFR has a scatter ot over an order of magnitude larger 
than what we observe at z~l. Given that the angu- 
lar momentu m properties are s hown to be similar to 
those at z <l (|Bouche et al.|2007| , if the sources targeted 
et al.l (|2009|) were less ( 



Forster Schreiber ct al. 2006 Gcnzcl 



et al.||2008 ), then a tight iW*-TFR would be delayed un 



til disk properties are homogenized across the observed 
populations — similarly for a transition from dynami- 
cally hot objects to disks which are primarily rotationally 
supported. Any inhomogeneity, including differences or 
changes in: 

• rotational versus dispersion support, 

• baryonic mass-to-light ratio across the disk, 

• baryonic versus dark matter fraction, 

• specific star formation rate, and 

• gas fraction or gas reservoir size and distribution 

across redshift bins should lead to an evolution in the 
M*-TFR. While various combinations of the above char- 
acteristics are likely correlated, they may not evolve in 



by Cresci et al.| Q2009p were less efficient star-formers 
than those studied here, perhaps this could explain the 
marked stellar mass deficiency indicated in Fig [7] How- 
ever, the earlier syste ms display higher ga s fractions and 
star- format ion rates ( Tacconi et al.|[2010 |, and consider- 
ing that they may have additional dynamical mass sup- 
ported by higher velocity dispersions only makes this dis- 
crepancy worse. Studies of less massive sources at z ~ 2 
and beyond as well as better constraints on the physical 
sizes of these components will be valuable in resolving 
t his possible d i screpa ncy. 

' 120111) note the 



Cresci et al. 



( 2009 ) and Gnerucci et al. 
decreased rotation-to-dispersion velocity fractio n (V/a) 
support in galaxies at z~2, as shown also in |F6rster| 



UUppUl U 111 tij <_4il ClJYl V I J CliU As ^5 IUJ U11U VV11 U'l..l w in 1 1 V/l u uvi I 

Schreiber et al. ( 2009[ ), as key in explaining the observed 
-/W*-TFR evolution. However this picture, which posits 
a greater degree of turbulent support in thicker disks at 
high redshift, should result in less rotational motion per 
unit stellar mass at earlier times, in contrast to the trend 
shown in Fig. 7J Given the negative zero-point evo- 
lution observedTa reduced V/a at higher redshift would 
thus imply a substantial evolution in the gas fraction. 
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Simultaneous changes in the dark matter fractions, bary- 
onic mass-to-light ratios, and disk geometries may lead 
to a less severe evolution in gas fraction. Understanding 
how these characteristics could conspire to keep the M*- 
TFR from more dramatic changes can be addressed by 
more sophisticated tools. 

Comparing the relative evolution in our study to that 
of semi-analytical models can thus provide interesting 
interpretations as to the predicted interplay between 
baryons and dark matter ove r these epochs. We note that 
the model of |Dutton et al.| ( 2011 ) departs significantly 
from the observed trends before z~l, whereas the Benson 
(2012) Galacticus models with several more forms of 
baryonic feedback are a closer match. This suggests that 
understanding the subtleties in the baryonic physics of 
galaxies is vitally important to our understanding galaxy 
assembly via tools l ike the Tully-Fisher relation. In fact, 
Dutton et al. ( 2011 ) see much less evolution in their bary- 
onic TuUy-l 1 isher relation over the same redshift range, 
signalling the effect of evolving gas fractions in conjunc- 
tion with their stellar components. Since the evolution 
of the dark halo drives the self-similar evolution in disk 
circular velocity in their model, the evolution present in 
both their baryonic and M^-TF Rs is being driven by dar k 



(2008) 



matter, similar to the results of So merville et al. 
However, we find in |Miller et al.| ( |201l[ ) thaT baryons 
may be the dominant driver of dynamics within T2.2, and 
furthermore, the dominance of dark matter would stabi- 
lize the di sk against the inten se star formation which is 
observed ( |Genzel et al.||2008| . Because the role of gas 
in disk galaxies throughout their formation is not ad- 
equately understood across a wide mass range at high 
redshift, it is likely that further work and exploration of 
the evolution of the baryonic Tully-Fisher relation will 
help to explain why we see so little evolution in the M*- 
TFR since z ~ 1.7, especia lly since we suppose that gas 
fractions evolve with time ( Tacconi et al.| 2010|, as well 
as the dark halo response to baryonic matter. 
Our main findings can be summarized as follows: 

1. Using LRIS on Keck I and techniques discussed in 
Paper I, we determine reliable rotational velocities 
using [O II] emission for 42 of 70 targeted star- 
forming galaxies selected from the redshift range 
1.0 < z < 1.7. We test the reliability of our 
extracted velocities with reference to our data at 
higher resolution and lower redshift. 

2. We find that the stellar mass Tully Fisher relation 
for this subsample is well-determined with up to 
60% increase in s catter compared to th e local re- 
lation defined by jPizagno et ah] (2005). Consid- 
ering also the relations found in Paper I and lo- 
cally, we find little zero-point shift corresponding 
to AM* = 0.02 ± 0.02 dcx from z = 1.7 to z = 0. 

3. Concerning the 28 galaxies which were not included 
in the above analysis, 6 are unresolved both in 
their broad-band imaging and in their spectra, and 
14 are spectrally passive. Although the remain- 
ing 8 show extended, resolved broad-band imag- 
ing, the emission is insufficiently extended to re- 
solve by LRIS with the extended exposure times 
of this study. We consider the impact of excluding 



these sources on our derived TF relation and note 
their properties. 

4. The modest evolution seen in the TF relation over 
< z < 1.7 contrasts markedly with results emerg- 
ing at z > 2. At face value it seems there is a 
dramatic change in the kinematic properties over 
a very short period of cosmic history however a 
wider range of sample stellar masses is required at 
z > 2 to understand the validity and significance 
of this result. 

While this paper was under review for publication, 
a study appeared on arXiv.o rg [astro-ph] pre-pri nt ser- 
vice, accepted for publication ( Vergani et al. 2012 1, which 
plots 23 galaxies from the MASS1V survey onto a stellar 
mass TF relation with similar or larger scatter as that 
found in our relation (cr int =0.32 or 0.52 and rms=0.48 or 
0.72) and little to some evolution in the stellar mass zero- 
poi nt (-0.05 ± 0.16 dex o r - 0.36 ± 0.11), depend ing on 



the |Bell fc de Jong| ( |2001 | or |Pizagno et al. (|2007|) slope 
adopted, respectively While the kinematics data was 
obtained on the SINFONI IFU, only 5 of these objects 
were observed using AO (and only 2 in the sub-sample 
used to plot the stellar mass TF relation). In their TF 
constraints cited above, they exclude over half of their 
sample where significant emission was detected, not only 
when the emission is unresolved, as in our study, but 
also when their velocity field is not well described by a 
symmetrically rotating disk, the kinematic position an- 
gle differs significantly from the morphologically derived 
major axis, or when their measured v/a < 1 (all crite- 
ria which are not included in our study). The addition 
of the modeled velocities in these excluded galaxies adds 
consid erably to the apparent scatter in the lVergani et al.] 
(2012 1 stellar mass Tully-Fisher relation, whereas our in- 
clusion of merger-like systems or galaxies where rotation 
curve centers appear offset to the HST-derived centers 
do not appear to result in the same large scatter. This 
may be due to their use of V m ax, described in this work, 



Paper I, and e.g. Courteau (1997) as potentially lead- 
ing to the underestimation of the true intrinsic maximal 
rotation velocity because of its dependancy on the distri- 
bution of emitting gas in the disk, inconsistent between 
galaxies in a diverse sample such as theirs. Also we find 
that assuming an intrinsic emission distribution based 
on broad-band optical or NIR HST imaging can often 
lead to systematic offsets in \ 2 minimization dynamical 
modeling, amplified by the scale of seeing and dispersion 
in the spectral data, as discussed in Paper I of this se- 
ries. If this effect is significant for slit-based modeling, 
then it would be just as damaging if not more so in 2-D 
IFU-based modeling. 

If the increased scatter we find in the TF relation in 
this study from 1 < z < 1.7 is truly intrinsic, then we 
should be able to determine what is different in those 
systems which have not yet arrived on the TF rela- 
tion. In future work, we will more closely examine the 
various characteristics of disk galaxies in our samples 
from 0.2 < z < 1.7: their estimated mass components, 
their dynamics, and their bulge-to-disk morphology. The 
full release of the WFC3 data in the CANDELS fields 



(Koekemoer et al. 2011) will allow us to compare our 



current bulge-to-disk decompositions to those in rest- 
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frame optical at our highest redshi fts observed. For ex- 
ample, it has bee n postulated (i.e., Bournaud et al.|2007 



Jones et al. 20101) that the development of a central bulge 



through inspiralmg, unstable clumps is key to stabilizing 
the disk. Alternatively, morphology could be driven by 
halo accretion type (hot versus cold), in which the cir- 
cular velocity shows little dependency o n morphology in 
the recent work by |Sales et al.| ( |2011 ). By examining 
the M*-TF relation as a function of the bulge to disk 
ratio and morphology, such hypotheses can be directly 
tested. Furthermore, an in depth comparison between 
semi-analytic and hydrodynamic models could provide 
valuable insight on the observed trends. In doing so, we 
aim to determine what is driving the emergence of the 
Tully-Fisher relation we observe by z ~ 1.7, and poten- 
tially, what in the formation processes of L* disk galaxies 
are leading them to largely stabilize by z ~ 1. 
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APPENDIX 
FURTHER SAMPLE DETAILS 



Here we more carefully detail the content of our Tully- 
Fisher sample, suggesting that we have a representative 
sample of rotationally supported galaxies at 1 < z < 1.7, 
including irregular and peculiar/interacting galaxies. 

For our initial selection based on visual morphology, 
the purpose for excluding unresolved, compact sources is 
that their dynamical measurements from emission lines 
are similarly unresolved (and in fact often non-existent). 
This accounts for an effective removal of about 1/5 of ob- 
jects, and an additional de- weighting of 1/10 of objects 
which meet the magnitude cut. As seen in the morpho- 
logical montage in Fig. [HJ the TF sample includes 4 
objects (1/10 of sample) that would be morphologically 
classified as early types or nearly unresolved (compact), 
as well as at least 8 objects (1/5 of sample) which would 
be considered peculiar, interacting, or undergoing a mi- 
nor merger. This fraction is consiste nt with the evolving 
merger rate fraction at this redshift ( Cassata et al.|2005 



Conselice et al.|2008| [Bridge et al.|2uTo| 



whether or not 

the initial visual morphological selection was biased to- 
ward mergers 



A d i fferent picture is presented in wo rks by Flores 
ci "all ( |2006| ); Kassin et al] ( |2007| ), and IPuech et al. 
(|2U08 ) which cite or suggest higher interacting and merg- 



ing fractions (at least 45% of their samples with com 
plex kinematics, with an additional 28% considered to 
be perturbed rotators in the latter work), which is used 
to explain increased scatter in the Tully-Fisher relation 
in these works. While it is well known in local studies 
that peculiar and interacting galaxies add some scatter 
to the idealized Tully-Fisher relation of a homogenous 



population (e.g., Tutui fc Sofue 



2002 Mendes de Uiiveira et al 



1997 Kannappan et al. 
W03p , this level ot scat- 



ter is sim i lar to what we find in our work. Rather |Flo 



res et al.| ( |2006[ ); |Kassin et al.| ( |2007[ ), and |Puech et al. 
] 2008 ) argue for the important contribution of peculiar 



and interacting galaxies on the larger scatter in their rela- 
tions. Although methods for identifying mergers vary in 
more complete sam p les at the same red shifts, e.g., |Kar^] 
taltepe et al.l (|2007|); iJogee et all 120091) (it is not trivial 



to compar e these due to d ifferent timescales of merger 
signatures, Conselice 2009), we nevertheless believe that 
measurement uncertainty from decreased signal-to-noise 
and resolution contributes to increased scatter in pre- 
vious intermediate-z studies, rather than an increased 
scatter arising primar ily from disturbed or complex kine- 
matics (see also, e.g., |Kannappan fc B arton 2004 Miller 



et al. 2011) 



Miller et a 



We h ave attempted in this work and in 
(12011) to clarify this important difference 



with extended integrations on spectra and a well-tested 
dynamical modeling method for deriving intrinsic rota- 
tion velocities. 

The three o bjects of overlap be twee n our Paper I sam- 
ple and that of |Flores et aT| ( |2006| ) and |Puech et al | ( |2008 1 
are all in the m ost kinematically complex classification 
of the Flores et al. ( 2006 1 scheme, so this suggests that 
we are not removing objects which may contribute max- 
i mally to the scatter seen in the TF relation, according 



to 



Puech et al. 



2008). Instead, as mentioned before, we 
see a similar amount of scatter as local studies which are 
inclusive to irreg ular, perturbed, and peculiar galaxies, 



Reyes et al. (i.e., 2011 ), which has a sample constructed 
to be representative based on SDSS. 

We find in Paper 1 for 0.2 < z < 1.3 galaxies that 
with only 1 hour of data, 1 /6 of the sample would have 
been in our spectrally compact category and 1/3 of the 
sample would have been considered passive in emission, 
i.e., not detected at all by DEIMOS. So over half of 
the sample that we include in our Tully-Fisher relations 
would be prematurely excluded in samples that make a 
significant emission cut based on ~1 hour redshift sur- 
veys from similarly sensitive instruments/telescopes as 
DEIMOS/Keck, such as the [Oil] emissio n criteria used 
to construct the MASSIV survey sample (Contini et al. 
201l | lyergani et al.|[2012] ). 

While in Paper 1 we implement a slitlet position angle 
offset (from major axis) velocity correction, significant in 
only 10 (8%) of objects, no such corrections were found 
necessary in this study (see Fig. [8b. 



z=1.21 
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Fig. 8. — Images of the galaxies in our TF sample, with either the HST F814W or F850LP and their 1 arcsecond slits overlaid. We also 
include the redshift of each target. 



